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Abstract 

A new model for elucidating the mathematical foundation of plasticity yield criteria is pro- 
posed. The proposed ansatz uses differential geometry and group theory concepts in addition 
to elementary hypotheses based on well-established experimental evidence. Its theoretical de- 
velopment involves the analysis of tensor functions and provides a series expansion which allows 
the functional stress-dependence of plasticity yield criteria to be predicted. The theoretical 
framework for the model includes a series of spatial coefficients that provide a more flexible 
theory for in-depth examination of symmetry and anisotropy in compact solid materials. It 
describes the classical yield criteria (like those of Tresca, Von Mises, Hosford, Hill, etc) and ac- 
curately describes the anomalous behaviour of metals such as aluminium, which was elucidated 
by Hill (1979). Further, absolutely new instances of stress-dependence are predicted; this makes 
it highly useful for fitting experimental data with a view to studying the phenomena behind 
plasticity. 

Keywords: geometrical model, analytic functions, manifold, anisotropic material, elastic-plastic 
material 

1 Introduction 

An accurate description of the structure, formation and behaviour of a solid clastic-plastic material 
requires the knowledge, among other facts, of the limiting stress it can withstand before it becomes 
plastic [32, 33, 15, 18, 2, 22]. In fact, plasticity concepts are widely used in a number of scientific 
and engineering field applications (in materials science, physics of solids, mechanical engineering, 
aeronautical engineering, geophysics, biomechanics and chemistry, among others) [4, 8, 1]. 

Plasticity involves a series of irreversible, history-dependent processes by effect of which a material 
develops fluency at a micro-, meso- or macroscopic scale in its transition from an elastic behaviour 
to a plastic behaviour [26]. Plastic processes involve plastic dissipation (i.e., the irreversible release 
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of stress or energy with energy transfer in the material). In fact, plasticity involves a variety of 
processes at different spatial scales having an also different associated grain size. In addition to its 
physical reality to plastic processes, grain size defines characteristic spatial scales for macroscopic 
plasticity in materials [19]. 

Unlike liquids and gases, solids are highly ordered systems, contain a vast amount of internal 
information and exhibit a high correlation among its constituent elements. As a result, plastic 
microscopic processes in solids (movement of dislocations, defect, etc) are usually relatively complex 
[19]. However, such processes exhibit some macroscopic symmetry by virtue of the solids structure 
and the physical laws they obey [4, 26, 5, 6]. This facilitates the macroscopic examination of solids 
by using a combination of differential geometry and group theory [3, 30]. Thus, a plastic process 
can be interpreted as a series of local transformations that possess some symmetry and provide local 
information useful with a view to establishing a global statistics for a solid material. 

The aim of this paper is to develop a unified approach of plasticity yield criteria that uses 
elementary hypotheses based on well-established experimental evidence. This study is based on the 
analysis of the transformation properties of the Cauchy stress tensor under orthogonal mappings, 
where arguments from the theory of Lie groups are applied. Also, it discusses the physical arguments 
for application of each criterion in relation to specific properties of the material concerned. This can 
be useful with a view to developing new criteria to address some special mechanical properties of 
materials (like anisotropy, hardening-softcning, etc). 

This paper is organized as follows. Section 2 describes some topics about plasticity in the 
framework of the internal variables theory. Section 3 presents a short description of classical plasticity 
criteria. In section 4 we show the postulates and approximations of the unified theory. From section 
5 onwards we develop this unified theory. Finally, a set of conclusions is shown. 

2 Plasticity function and plastic potential. Flow rule 

Within the framework of the Internal Variables Theory (or hidden variables), an inelastic solid is 
one in which the strain at any point of the solid is completely determined by the current stress 
and temperature there plus a set of internal variables [26, 28]. The internal variables (scalars or 
tensors) have physical or mathematical meaning and allow to complete the internal description at 
any point of the solid (for example: the past history of the stress and temperature at the point, 
large-deformation plasticity, hardening and softening, structural and induced anisotropy, etc) [26, 
29] . Thus e is a function of the material state (a, T, £) at any point 

e = e(a,T,0 (1) 

where e denotes strain variables, a stress variables, T the temperature, and £ some internal variables. 
Additionally, the rate of evolution of the internal variables £ is determined by the state 

i = fa,T,0 (2) 

known like equations of evolution or rate equation for internal variables. 

For inelastic solids it is generally assumed that strain variables can be decomposed additively 
into elastic strain e e and inelastic strain e l [8] 
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e = e e + e 1 - (3) 

where the inelastic strain occurring in rate-independent plasticity is usually denoted by e^rather 
than by e\ and is called the plastic strain. 

In the context of internal variables theory and rate-independent plasticity, a plasticity yield 
criterion consists of a series of mathematical conditions mutually relating stress, temperature and 
internal variables, which define the material states (a, T, £), where the point of the solid concerned 
becomes critically plastic. A critical state in this context is a state where the elastic-plastic material 
starts to yield plasticity. Therefore, these conditions constitute a boundary (plastic limit) between 
the elastic and plastic state in the solid material (at points). Mathematically, plasticity yield criteria 
are formulated in the following general form: 

/(<r,T,0=O (4) 

where / denotes the plasticity function at a given plastic state (a, T, £). Therefore, the plasticity 
yield criterion of a solid (at points) defines the stress multiaxial states where it will yield critical 
plasticity; the set of such states describes the criterion surface of the solid f4} . 

Plasticity can be defined in basic terms by using some approximations that usually hold in 
practice. Thus, the influence of the temperature T on solids at a constant, ambient level is usually 
negligible provided they are scarcely sensitive to changes in this variable and far from their melting 
point (temperature-independent plasticity). The influence of strain is usually negligible if we consider 
all viscoplastic processes to be "infinitely" slow compared to the material relaxation time (rate- 
independent plasticity) [26] . Time- independent plasticity needs usually to be considered if the target 
is not the time-evolution of the solid. In the plastic limit, an index P is used to denote parameters 
and quantity values in such a limit. The basic plastic unit in a solid is the macroscopic spatial point, 
which can be equated to the physical concept of grain but need not coincide with it. Also, because 
a solid usually exhibits high correlation among its elements, characterizing each point in it requires 
defining its correlation with its neighbourhood. As in the theory of elasticity, this entails describing 
the stress state at each point in terms of a second-order symmetric tensor (2-tensor) a = <xy , called 
the Cauchy tensor, the symmetry of which arises from the stress equilibrium relation at the point 
in question [24]. Therefore, in these conditions, the plasticity yield criterion ((4]) at each point in a 
solid can be defined as follows: 

/(^■;0 = o (5) 

which represents the criterion surface in the stress space. 

The plastic potential, is a measure of smoothness (differentiability) and convexity of /. Thus, 
if the plasticity function / is smooth (differentiable) and convex in stress space, then it will coincide 
with a specific plastic potential (/ = $). The plasticity function / = $ is convex if for a given 
stress and strain rate condition cry and ifj, any stress a*j inside or on the criterion surface obeys 
the following relationship [26] : 

(a t] - a* 3 )e[ 3 > (6) 

where efj denotes the plastic strain rate. Specifically, if / = <f> is twice differentiable, then $ will be 
convex if, and only if, its Hessian matrix, defined as 
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H a = nTiir" (7) 



<9 2 $ 

da id a j 

is positive scmi-dcfinitc (i.e. if the eigenvalues of Hij arc only positive or zero). When the eigenvalues 
are only positive non-zero (positive definite matrix), / = <i> is strictly convex. We have used here 
— as we have throughout — the Einstein summation convention for repeated indices. 

If / is smooth (differentiable) and convex, then ef^ will be unique in any plastic stress state 
crij . Under these conditions, the function / = <f> can be assigned a flow rule that gives the plastic 
strain rate. Within the framework of the plastic potential theory of Von Mises (1928), each plastic 
potential $ has a flow rule that is associated with the plastic potential. This flow rule is given by 

where <7y is the Cauchy tensor, ef^ the plastic strain rate tensor or plastic flow tensor and A a positive 
plastic multiplying factor. Based on eq. ©, if / = <f>, then the plastic strain rate efj at the very 
start of plasticity will be normal to the criterion surface / = <&; this constitutes the so-called plastic 
flow normality rule. Hecker [13] conducted a systematic study of a large amount of experimental 
data in metals and found that the normality rule was never broken. However, there is evidence that 
flow normality rule does not hold in soils, granular materials, etc (materials with non-associated flow 
rules, where df/daij is not proportional to d^/daij) [8, 26]. 

Defining rate- independent plasticity requires introducing the concept of plastic dissipation, which 
is expressed as 

d = rr, , (9) 

The parameter d is a measure of the power per volume that is lost (dissipated), usually as heat, 
through deformation. The plastic dissipation as defined in eq. © for the flow rule associated to / 
will be maximal if, and only if, the function / = $ is convex this is the so-called principle of 
maximum plastic dissipation [31, 16]. 

Establishing the total dissipation over a given time interval entails defining in terms of time- 
dependent plasticity the amount of irreversible work per unit volume due to stress as follows: 

X (t) = W p (t)= f a l3 (t)e?(t)dt (10) 
Jo 

x(t) is also known as the work-hardening parameter (an internal variable), which is a scalar quantity 
and dependent on the plastic time-evolution of the particular solid. 

For the study of the hardening and softening of the solids materials is necessary the evolution of 
plasticity function with the internal variables, £ = £ Q , given by [26] 



/ 



a— con 



St., T=const. ' o£ a O^a 



where H > for hardening materials, H < for softening materials, and H = for perfectly plastic 
materials (in the latter case / is independent of the £ a ). Here A is a positive continuous function of 
state variables. 
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Based on the foregoing, for the theory to be properly addressed, the plasticity function / should 
be convex in the stress- related variables (/ = <&). The dependency of / on £ describes the anisotropy 
(structural and induced) of materials, while the dependency of / on £ describes the hardening- 
softening of materials. 

3 Plasticity yield criteria 

The following section describes the most relevant plasticity yield criteria, with emphasis on their 
particularities (see [32, 33, 15, 18, 2, 22]). 

The Tresca criterion [32] which is among the earliest plasticity yield criteria, can be expressed 
as a unique function of the algebraic invariants (Ji, J2, J3) of the stress deviation tensor J. Based 
on it, a solid will become plastic when it reaches a multiaxial state where the multiaxial tangent 
stress equals the critical uniaxial tangent stress, a p . This criterion can be expressed as a completely 
diffcrentiable relation 

fi&ijl&p) = [{o-xx - o yy f + Aa 2 y - Aa 2 P ] [(a vy - o zz ) 2 + 4a 2 yz - Aa 2 P ] 

x[(<J zz - o xx ) 2 + 4a 2 xz - 4a 2 P ] (12) 

based on which plasticity can only be reached on independent planes, as revealed by factoring the 
total (volume) multiaxial state into its partial (surface) multiaxial state. In other words, plasticity 
develops on planes. This criterion exhibits good agreement with experimental results for certain 
ductile metals. At microscopic level the movement of dislocations along slip planes is responsible 
for permanent deformation [4]. Note that a P depends on the internal variables, £, i.e. a P = cr P (£). 

The Huber-Von Mises criterion [33, 21] is among the most widely rule used in this context. 
It is also known as the J2- criterion since it is formulated as a unique function of the algebraic 
invariant J 2 of the stress deviation tensor. Based on existing experimental evidence, this criterion 
is applicable to ductile materials. Hencky [14] provided an energy-based interpretation by which 
the critical plasticity state is reached when the distortion energy per volume (i.e., the deformation 
energy per volume in the absence of volume changes) in the multiaxial tangent stress state equals the 
distortion energy per volume in the critical uniaxial normal stress state. This criterion is isotropic 
and can be formulated as follows: 

f(0-ij-,0-p) = (0 XX - Oyy) 2 + {0-yy ~ O Z Z f + (0 ZZ ~ (T XX ) 2 + 6 (c^ + O yz + (7^ ) ~ 2<7p (13) 

Stress-wise, the Von Mises criterion indicates that critical plasticity is reached when the modulus 
of the multiaxial tangent stress equals that of the critical uniaxial normal stress, a p . Note that a P 
depend on the internal variables, £, i.e. a p = o p {^). Unlike the Tresca criterion, the Von Mises 
criterion does not factor the plasticity function into multiaxial plane stress components; rather, it 
assumes a mutual dependence among the multiaxial stresses, which leads to a volume plasticity 
expression. In ductile materials, which meet the Von Mises criterion quite well, plasticity results in 
distortion by effect of flow processes occurring with virtually no volume change. 

The Hill's anisotropic criterion or Hill's first criterion [15, 17] provides a general description of 
materials with anisotropy (whether structural or induced) and orthotropic symmetry (i.e., materials 
where each point possess three mutually normal planes). Because each plane in an orthotropic 
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material can be defined in terms of only two parameters, Hill's criterion can be formulated in terms 
of six independent parameters. This constitutes a generalization of the Von Mises criterion to 
anisotropic materials of the form 

A(a vv - <j zz f + B{a zz - <j xx f + C(a xx - oyyf + 2Da 2 yz + 2Eo 2 xz + 2Fa 2 xy = 1 (14) 

where A = A(£),B = B(£),C = C(£),D = D(£),E = E(£),F = F(£), depending on the internal 
variables, are constants defining the degree of anisotropy in each direction and can be expressed in 
terms of the Lankford coefficient. In addition, the directions x, y, z should be the principal anisotropy 
directions for the material — otherwise, the Cauchy tensor should be transformed as required in order 
to have it coincide with the principal directions. In fact, eq. ([M]) is a particular case [26] of 

&zjB i jkiO'ki = Op (15) 

which is the more general quadratic form for the components of the stress tensor as defined in terms 
of the 4-tensor Bijki. This tensor fulfils the symmetry conditions for 4-tensors in the theory of 
elasticity, which are given by 

Bijki = Bfruj 
Bijki = Bjif-i 

where the first relation defines the symmetry by joint exchange of index pairs ij «-> kl and the 
second the symmetry by exchange between index pairs of the type i <-> j and/or k <-> I. This 
criterion considers no effects of the mean stress (hydrostatic pressure) on plasticity; therefore, one 
must introduce the additional condition Bijkk = 0, which allows the quadratic form of such effects 
to be neglected. 

The Hosford criterion [20] and Logan-Hosford criterion [25] are two generalizations of the Von 
Mises criterion that ensure convexity by introducing a parameter m > 1. The Hosford criterion, 
which is applicable to isotropic materials, is defined as 

*p) = \\°i- ^ + \ ka - <ra\ m + \Wx~ ^1™ - (17) 

where oi, 02, 0-3 are the principal stresses of the Cauchy tensor. Note that a P depends on the internal 
variables, a P = cr P (£). The Logan-Hosford criterion is a generalization of Hill's anisotropic criterion 
(|14p that is used to describe anisotropic materials and expressed as 

A' lax - a 2 \ rn + B' \* 2 - o- 3 \ m + C ki - o- 3 \ m = 1 (18) 

where the constants A' = A'(^),B' — B'(£),C = C"(£), depending on the internal variables, 
constitute measures of anisotropy in each direction and can be expressed as a function of the Lankford 
coefficient. 

Hill's generalized anisotropic criterion or Hill's second criterion [18] provides an accurate descrip- 
tion of the anomalous behaviour of some metals such as aluminium [37]. The criterion is expressed 
in terms of principal stresses of the Cauchy tensor: 

a \<r 2 - cr 3 1 m + b cr 3 — cri| m + c |cti - o- 2 \ m + d \2o\ —01 — 03!™ + 

e \2a 2 - ai -a 3 \ m + f \2<j 3 -a x - a 2 \ m = (19) 



(i 



where m G Jf is a parameter that must fulfil the condition m > 1 for the criterion surface to be 
convex. Also the constants a = a(£),b = &(£),c = c(£),d = d(£),e = e(£),/ = /(£), o> = cr p (£) 
depending on the internal variables, £, and which can be expressed as a function of the Lankford 
coefficient, measure anisotropy in each principal direction. For example, if the principal directions 
1-2 in a material define a symmetry plane with mutual isotropy and are anisotropic with respect to 
direction 3, then a = b and d = e (planar isotropy). On the other hand, if all three directions are 
isotropic, then a = b = c and d = e = f (complete isotropy). 

In the plastic potential theory of Von Mises [34], each plastic potential $ is assigned a flow rule. 
The plastic potential for the particular case of the Von Mises criterion is J2 and its associated flow 
rule defined by the Levy-Mises equations. Koiter [23] developed a generalization of the previous 
theory where the plasticity function is defined by a series of plastic potentials $j each having an 
associated flow rule of the type described by eq. © above. If the plasticity function / is expressed 
as a linear combination with positive coefficients of the potential functions $ j (convex) , then / will 
be convex. 

The plastic anisotropy description of Barlat et al. [2] originated from an isotropic plasticity 
function. Structural anisotropy in the material was introduced via a series of linear transformations 
represented by a 4-tensor acting on the Cauchy tensor, the latter itself acting on the anisotropic 
material. By effect of the transformations, the stress tensor absorbs structural anisotropy of the 
material. 

The Karafillis-Boyce criterion [22] uses a convex combination of two plastic potentials as plas- 
ticity function. The potentials are based on Hosford's isotropic criterion (fl7|) and their degree of 
mixing is adjusted via parameter c 6 [0, 1]. Thus, the isotropic criterion is formulated as 

om 

(1 - c) [\ Sl s 2 \ m + \s 2 - s 3 \ m + 1*3 - Sl H + [\si\ m + \s 2 \ m + = 2a£ (20) 

where m is positive and non-zero, and si, S2, S3 are the principal values (eigenvalues) of the stress 
deviation tensor. Note that a P , c depend on the internal variables, o~ P — cr p (£), c — c(£). This 
criterion is a particular case of Hill's second criterion (fT9f . The plasticity yield criterion ([20|) repre- 
sents an isotropic, convex criterion that can be made anisotropic by applying linear transformations 
representing a 4-tensor acting on the Cauchy 2-tensor, which in turn act on the anisotropic mate- 
rial (Barlat et al., 1991). Introducing appropriate symmetries of a material in the transformation 
4-tensor allows all possible states of structural anisotropy in the material to be considered. The 
Barlat and Karafillis-Boyce criteria rely on the theory of representation of tensor functions [35] . 

Those criteria that consider hydrostatic pressure dependence assume plasticity in some materials 
including metallic foams and polymers to be a function of the hydrostatic pressure o~m acting on 
the solid. This effect has been considered by using various general criteria such as those of Drucker- 
Prager [10, 11], Caddell [7] and Deshpande [9], which are modifications of the Von Mises and Hill 
criteria including a certain dependence on the hydrostatic pressure (o~m)- 

The plasticity yield criterion is established from the set of microscopic and macroscopic properties 
of the material. As a result, the formulation of each criterion depends on a combination of parameters 
of the material describing its anisotropy, crystal structure and hydrostatic pressure-dependence, 
among other properties. There are three general types of models for analysing plasticity, namely: 
microstructurc or macrostructure and mixed. Micro structure models establish plasticity yield criteria 
from the microstructure of each material. On the other hand, macrostructure models, also referred 
to as phenomenological models, rely on a phenomenological analysis of the macroscopic behaviour 
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of each material to establish such criteria. Finally, mixed models are combinations of the previous 
two and usually provide the more accurate descriptions of plasticity. 

4 A unified theory of plasticity 

The objective of this paper is to develop a new macroscopic theory that unifies the plasticity yield 
criteria. This theory is based on postulates well-established from experimental data and theoretical 
considerations. The used method is based on orthogonal Lie Groups to describe the classical isotropic 
yield criteria; an increase of the symmetry group allows to consider classical and new anisotropy 
yield criteria in the solid materials (new mechanical properties). 

At microscopic level, the plasticity in a solid is produced by a set of slips, generated by the 
movement of dislocations, defects, etc [19]. The stress acting on the solid is the generator of these 
slips. We want to take into account all these microscopic dynamical processes at the macroscopic 
scale. So the theory is based on the decomposition of a solid material into parts as macroscopic 
points p (Figure 1). These macroscopic points do not need to have any physical reality (grains 
or others), but points have macroscopic information about microscopic dynamical processes. This 
macroscopic information is purely statistical, so a macroscopic point is a statistical concept that 
takes into account all the microscopic information inside it. The set of slips inside a macroscopic 
point induces the irreversible macroscopic movement of this point (with a well determined spatial 
directionality). The global movements of the points result in macroscopic plasticity that follow 
macroscopic laws (theory postulates). So, the proposed unified theory examines plasticity at each 
macroscopic point p in a solid. The plasticity of the macroscopic points are described by orthogonal 
transformation groups (directionality of plasticity) acting on Cauchy tensor (generator of plasticity) 
in a convex manner (measure of plasticity) . 

The macroscopic anisotropy description by coordinate tensors appears in the infinitesimal or- 
thogonal transformation of Cauchy stress tensor (which acts on the solid point) with respect to 
a reference system (for the solid point). The directions of these infinitesimal transformations are 
limited only by the internal structure (microscopic dynamic processes) of macroscopic solid points. 
These directions of transformation (anisotropy characterization) are described by the coordinate 
tensors (parameters that depend on the internal variables). A measure (mathematical norm) of the 
infinitesimal orthogonal transformation of Cauchy tensor gives the macroscopic effect of plasticity 
(plasticity function). The ensuing point information can be used to obtain a global points statistic 
for the whole solid. 

Certain approaches are considered in this theory. The potential effects of temperature and 
strain arc ignored (temperature- independent and rate- independent plasticity), and so is the time- 
dependence of all parameters and variables since the only target is critical plasticity in the solid 
(time- independent plasticity). Our theory relies on the following three postulates: 

1. First postulate: the macroscopic units or points that become plastic in a solid are described 
by second-order Cauchy tensors. 

2. Second postulate: hydrostatic-pressure independent plasticity acts at every point in the 
Tangent Space to the Cauchy Tensor (TSCT). Any type of plasticity at least needs a component 
in this space TSCT. 

3. Third postulate: the plasticity function is convex in the stress space. 
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Within the framework of the proposed theory, solids are represented in a three-dimensional Euclidean 
space -E© where each point or neighbourhood p G E© is defined by the pair p = {p x , <r kl } £ E$5§, 
with indices X,k,l = x, y, z. The set of points is referred to E^ via a coordinate system or reference 
system {O, x x = (x, y, z)} that is shared by the whole solid; the coordinate system has O as its origin 
and x x = (x,y, z) as its axes. Therefore, a point p in the solid lies at p x in the reference system 
and has the symmetric second-order stress tensor or Cauchy tensor acting on it, i.e., a = a kl {first 
postulate); this postulate is valid only for solids with "near-action" internal forces [24]. Since the 
solid is examined point-wise, each point p possesses its own reference system, {O p ,x x = (x, y, z)}, 
with its origin at p x and x x — (x, y, z) as axes. So the set of stresses acting on the points in the 
solid define a 2-tensor stress field in the space E§3§, which endows the solid with geometric structure 
(Figure 1). 

The reference system for each point, {O p ,x x = (x, y, z)}, is arbitrary only when the ensuing 
equations are invariant under a reference system change (tensor equations), that is, independent 
of the particular reference system. From the geometric point of view this guarantees that the 
tensor equations that describe the solid contain, codified as geometric information, all the physical 
behaviour throughout the directions of reference system changes for the corresponding space. 

If the Cauchy tensor is invariant under a change in the origin and axes over the solid, {0 Pl , x x = 
(x, y, z)} —* {0 P2 ,x x = (x,y,z)}, then the solid is stress-homogeneous. If the Cauchy tensor is 
invariant under a change in the axes directions, {O p ,x x — (x,y,z)} — » {O p ,x x = (x,y,z)}, then 
the stresses at the point p in the solid are isotropic. The degree of homogeneity in the solid can be 
assessed by examining it point-wise and using the results to develop a global statistics for the entire 
solid. In this work, we will assume the behaviour of the solid is similar to its individual points or, 
in other words, that examining a single point will be the same as studying the entire solid. 

The intrinsic properties of the solid at each point are examined via a series of transformations 
with a fixed origin on it. A transformation of the Cauchy tensor, a kl , at point p can be interpreted 
in two completely equivalent ways. In the first one, the reference system x A has the Cauchy tensor 
o~ kl ; therefore, a transformation x A — > x A in the reference system will lead to a new reference 
system x A with a Cauchy tensor cr* 3 '; this constitutes a reference system transformation or active 
transformation. In the other interpretation, the reference system x A remains unchanged and the 
Cauchy tensor transformation, a kl — > ct 2 - 7 , produces a new tensor that is expressed (oriented) in the 
reference system 5c A ; this is a Cauchy tensor transformation or passive transformation. These two 
interpretations are mutually related by a negative sign. 

The proposed theory relies on a continuous linear transformation of the stress at each point p 
represented by the tensor AV, ■ Such a transformation does not alter the origin of p, which is defined 
by the coordinates p A , but only its coordinate axes, x x . Thus, the transformation operator A 1 ^ 
acts on the stress tensor a kl to give the transformed tensor <r y in accordance with the following 
expression: 

A%a kl =av (21) 

where the multilinear transformation AS is a 4-tensor of indices i,j,k,l — x, y, z with two covariant 
(k,l) and two contravariant(i, j) indices. Mathematically, it follows that this expression is a tensor 
equation (i.e., that it is invariant under a coordinate system change). Since the resulting transformed 
tensor, , has a fixed origin (the same that o~ kl ), the transformation can be tangent or normal to 
the tensor a kl . In the former case, it will produce rotations and/or reflections of use for studying 
anisotropy; in the latter case, it will result in dilatation or dilatation-reflection useful with a view 
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to examining hydrostatic pressure-dependence. 

The plasticity at a point p under the influence of an arbitrary initial Cauchy tensor depends on 
the tangent and normal spaces to the Cauchy tensor. These spaces arise directly from the geometric 
generalization of the two types of components of the Cauchy tensor. The direct sum of the tangent 
and normal spaces constitutes the overall stress space at the point concerned. As implied by our 
second postulate, a solid will yield plasticity at a given point only if a component on the tangent 
space exists at such a point. As a result, plasticity starts in the associated pure tangent space 
(hydrostatic-pressure independent plasticity) or, in other words, it never starts in the associated 
pure normal (hydrostatic- pressure dependent plasticity) space. In fact, there is solid experimental 
evidence that plasticity never arises under conditions of pure hydrostatic pressure [4, 26]. Thus, 
some solids exhibit slight elastic (reversible) deformation rather than plasticity, even at very high 
pure compressive hydrostatic pressures. However, high stress-induced hydrostatic pressures in a solid 
under tension can lead to an unexpected spatial stress concentration eventually leading to fragile 
fracture [5, 6]. 



5 Rotational transformation of the Cauchy tensor 

Properly examining plasticity in a solid entails starting in the Tangent Space to the Cauchy Tensor 
(TSCT) associated to each point p in the solid, and applying a pure orthogonal transformation R^ 
such that 

R%a kl =a«, (22) 

where the orthogonality condition is described by means of the matrix equation R T R = I. Based 
on this tensor equation, the Cauchy tensor a 1 associated to each point p in the solid is rotated by 
R£, according to a given parameter, which provides the new associated tensor a lJ , which will be the 
new associated tensor acting on p. Stress rotations according to equation (|22[) provide a point-wise 
definition of hydrostatic-pressure independent plasticity in the solid. In this sense, rotations of the 
Cauchy tensor at p can be interpreted physically as due to some tangential stresses acting on the 
point. In this respect, the hydrostatic-pressure independent plasticity interaction (i.e., that having 
no hydrostatic component) is fully defined by the rotations W£, and their associated transformation 
parameters. 

The orthogonal rotation transformation R 1 ^ is a three-dimensional 4-tensor (i.e., one with 3 4 = 81 
components) equivalent to a Mohr transformation if parameterized in Euler angles Qe = 6e((», 0, 7) 
[27]. As noted earlier, the orthogonal rotation transformation is necessary for hydrostatic-pressure 
independent plasticity to develop at individual points in a solid. However, additional, non-orthogonal 
transformations (e.g., dilations) can also be applied that will alter the hydrostatic-pressure indepen- 
dent plasticity conditions imposed by a rotation (see Section 11). The remainder of this section, 
and all subsequent ones up to the eleventh, is devoted to examine hydrostatic-pressure independent 
plasticity as defined in eq. (|22|) . Section 11 is concerned with hydrostatic pressure-dependence. 

The transformation equation (f2"2")) connects geometric objects via double indices (i.e., 9 compo- 
nents). However, an equivalent representation connecting objects via single indices (i.e., 3 compo- 
nents) can be obtained by using the following relations: 

n /,>': uP^Rlu 1 (23) 
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u j = T ]q v q ; v? = T jq v q (24) 

where (|2"B1 defines the rotational transformations of contravariant vectors and (fM]) the general 
contravariant transformations for a vector under coordinate system changes (i.e., the definition of a 
1-tensor). From the rotational transformation (|23|) and general coordinate transformation ([24]) for 
each 1-tensor, one can obtain 

u l = //;//' = R)T^v q = RjjT^RPv. 

Since the coordinate system change occurs via the stress 2-tensor (Cauchy tensor) T- 79 = cr- 79 , we 
can express the transformation by means of a rotational operator 

fr« = R l k o kl R{ <la = R T oR, (26) 

the latter describing the transformation rule for the Cauchy tensor. It is straightforward to verify 
that this expression is equivalent to (|22|) . The 2-tensor rotations of the types R k and Rjc&n be 
adequately realized in terms of 3 x 3 matrices. Relation (|26[) is invariant under a general change of 
the coordinate system. 

Therefore, the transformations^ , R\ lead the plastic process at each point p in the solid. Their 
study is facilitated by considering the three-dimensional special orthogonal group .SO© of rotations 
in the three-dimensional Euclidean space E$3§. The main advantage of this approach is that SO (O 
possesses both the structure of a group of transformations and a difierentiable manifold, which 
enables us to confer additional invariance properties to the tensor operators defined over it [12]. 

The elements of the special orthogonal group .SO ([3]) are the orthogonal transformations of the 
Euclidean space E§3§ that describe orientation and length preserving movements. Such rotations are 
usually represented by orthogonal real 3x3 matrices R T R = I with unit determinant (det R = +1). 
and the group product operation is the usual matrix multiplication. Topologically, .SO® is a 
compact, non-simply connected group [30]. As Lie group, S'O® is simple, i.e., the only normal 
subgroups it contains are the trivial ones: itself and the identity group. In particular, its Lie algebra 
so([3]), which coincides with the tangent space at the identity element, is also a simple Lie algebra. 
As a consequence of the Lie structure the tangent bundle TSO§5§ = Up^eso© TSO^ pV inherits 
special properties that will be useful in our later analysis. 



I => f lp = R)T ]q R p q (25) 



6 Action of the group SO($) on the Cauchy tensor 

Using the adjoint representation of SO^, the action of the rotations R\\ in ((22]) or the equivalent 
pair R k , R\ in P5]l can be explicitly expressed in terms of the Euler angles 9^ = Qe(o, {3,^) as: 

a xx = <j xx cos 2 7 + a yy sin 2 ^ + 2a xy sinj cos 7 
a yv — a vy cos 2 7 + a xx sin 2 ^j — 2a xy siwy cos 7 
a zz = a zz 

< (27) 
a xy — a xy (cos 2 7 — sin 2 j) + (a vv — a xx )sinj cos 7 v ' 

a yz = a yz cos 7 — a xz sin'j 

a xz = a yz sin-f + a xz cos 7 
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a xx = a xx cos (3 + a zz sin (3 + 2<j xz sinf3cos (3 
a yy = a yy 

a zz — a zz cos 2 + cr xx sin 2 (3 — 2a xz sin(3 cos/? 
a xy = a yz sinf3 + a xy cos (3 
a yz = a yz cos/3 — a xy 'sin/3 



(28) 



'■(cos 2 (3 - sin 2 (3) + {a z 



")sin(3 cos/3 



a xx 
a yy 
u zz 

<7 X 



cos 2 a 



77 ■ ? 

t sin a 



2a yz sina cos a 



a zz cos 2 a + a vv sin 2 a + 2a yz sina cos a 
a xy cos a — a xz sina 



(29) 



a yz — a yz (cos 2 a — sin 2 a) + (a yy — a zz )sina cos a 



--xz _ xy 



a y since + a cos a 

A rotation at the point p therefore represents a transformation from the reference system 
{O p ,x A }, with a tensor <j kl , to the new reference system {O p ,x A }, with associated tensor a 1 ^ . The 
components of the Cauchy tensor in the rotated reference system cr y are expressed in terms of a kl 
by means of equations (|2T|) - (f2"9"]) . 

The rotation by angle Oe = @_e(ck, (3, 7) can be easily visualized by viewing the reference system 
in the direction from the positive ends of the axes to their origin and rotating x x anti-clockwise to 
the new system, x x , by an angle a, (3 and 7 about the x, y and z axis, respectively (see Figures 
2 and 3). Obviously, because of the non-Abelianity of SO§3§, the order in which the three partial 
rotations are performed strongly influences the general rotational outcome. 

The action of a generic rotation on the Cauchy tensor is the result of applying successive partial 
rotations (Figure 3) about each axis in a given sequence. Thus, applying the rotations following the 
sequence a — > (3 — > 7, we get 



hi 



(30) 



which is equivalent to the following expression for rotations in the group 50([3|): 

* ij = [ J R( 7 )i?(/3)i?(a)]l.a fci [R(a)R((3)R(l)]i * * = R T (l, P, a)aR(a, (3, 7) (31) 

with partial rotations performed following the sequence a — > (3 — > 7. We observe that, although 
equations (|27f - |29|) depend on the explicit choice of rotation parameters, the transformation does 
not alter the form of equations (|30|) and ([3~Tj) . which, as noted earlier, are invariant under a change 
of coordinate system. 



7 Infinitesimal action of SO ^ on the Cauchy tensor 

Because of the properties of Lie algebras and the local diffeormorphism on a neighbourhood of the 
identity element [12], the action (|3Tj) on the Cauchy tensor can be also described in terms of the 
infinitesimal rotations, which turn out to be the generators of soQ. 

The Lie group SOjl) and its associated Lie algebra, so®, are related via the exponential 
mapping. Given an element X v £ so((3]), the exponential map provides an element of SO ((3|) by 
means of the relation: 
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Exp(X^) = R(9 V ) — e . (32) 

Since SO§5§ is a compact group, it is ensured that Exp is a surjective mapping [30]. Differentiation 
with respect to the parameter and evaluation at the identity element provides the operator X: 

^ m*R(Q v )\ e „=o = Xv (33) 

d& R ( e E)\e=Q = P v M* R ( 9u )\e»=o = P VXv = x 

The operators X in the Lie algebra therefore correspond to infinitesimal rotations and can be 
expressed as linear combinations of the operators for the local basis X v = djd9 v in the Lie algebra 
so((3]). As a result, each infinitesimal rotation can be expressed as X — p v X v = p v (d/d9 l/ ), which 
has p v as its local coordinates. 

The actions of infinitesimal rotations on the stress field are determined by the exponential map 
(|3"2")l . which relates the group and the Lie algebra; this equation is a function of the infinitesimal 
generators X v of each partial rotation. A first-order series expansion about the angle 6 V — (group 
identity) yields the expression: 

R(0 V ) = e e " x » =1 + X v 50 v + ... (34) 

where Qe = 8s(a,/3, 7) = Qe(9 1 , 9 2 , 9 3 ) denote the Euler angles associated to the first-order 
infinitesimal rotations 89 v , which can take the values 89 v = 5a, 6(3, 5^. In addition, each infinitesimal 
angle 89 v has an associated infinitesimal generator X v = djdB v . Substituting (|34|) into the rotational 
expression for a 2-tensor [eq. (|3ip ] yields the expression describing the action of partial infinitesimal 
rotations, by an angle V , on the Cauchy tensor. For a first-order approximation, such an expression 
is of the form 

a(6 v + 56") = (I + X V 56») T a{8 v ) (I + X V 5B V ) = , . 

= (7(9") + {a[6 v )X v + XTa{9»)}59» + O 2 {86") 1 ' 

where eq. (|3T|) has been applied on the assumption that a — a(9 v + 89") and a — cr(6 )I/ ) (i.e., the 
neither the origin nor the rotation centre of the tensor is altered, since the transformation of the 
Cauchy tensor involves rotations only, so that the tensor remains unchanged but is evaluated at 
a different angle). Equation (I35|) can be used to establish the first-order variation at each partial 
angle, which is of the form 

a(9 v + 89") - a(6 v ) = a(9 v )X v 89 v + (X u f a{8")89 lJ & . . 

= aX v + {X v ) T a= [a,X u ] 1 ' 

where the infinitesimal generators X v for S'OQ are antisymmetric and the binary operation [a, X v ] = 
aX v — X v a is the commutator between the stress field and the infinitesimal generators X v . Equation 
(J3SJ) is a function of the specific basis X v associated to the Euler angles. Expressing it in terms of an 
invariant generic infinitesimal rotation X describing the tangent spaces to .SO© in any arbitrary 
coordinate system entails introducing the local coordinates p u associated to the local basis X u for 
the infinitesimal rotation X. 

This results in a first- order variation of the tensor stress field a with the invariant generator 
(infinitesimal rotation) X — p v X v of the form 
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da\ 1 = [a,X]=aX-Xa = p" — (37) 

which is the equation for a 2-tensor. Since the operator X is antisymmetric, then da\ 1 = [o~,X] = 
— [X, a] = — dXL, based on which the variation of a with X is equivalent to that of X with a with 
the opposite sign. This is consistent with the two interpretations of a transformation provided in 
Section 4 as applied to a first-order variation. 

The tangent space produced by the first-order variation (|37|) is a linear combination of Euler 
angles (or, similarly generators X u ) where each term is linearly dependent on the angles (or gen- 
erators); therefore, eq. (f3"T|) comprises a linear basis of generators X v (an Euclidean space). This 
linear approximation describes TSCT at a first-order local level only (i.e., in very small or local 
neighbourhoods) (see Figure 4). Expanding such neighbourhoods entails using higher-order approx- 
imations involving non- linear combinations in the X v generator basis (i.e., non-Euclidean spaces for 
the different orders) (Figure 5). As shown below by explicit calculations, all spaces approximating 
TSCT are traceless (i.e., the complete approximation are in fact tangent to the Cauchy tensor). 
This allows TSCT to be described by means of expansion series of different order. 

If one assumes anisotropy at each point in a solid to be more accurately described by including 
higher order-terms in the series, then the geometric interpretation of such anisotropy suggests that 
the more local the space — and the smaller the neighbourhood — , the higher will be the degree of 
isotropy. The inclusion of terms of an increasingly higher order gradually makes the local structure 
complex enough to facilitate extension of each infinitesimal neighbourhood, with improves approx- 
imation, to the entire TSCT. The complex behaviour of solids as concerning symmetry, anisotropy 
and various other properties dictates the structure of TSCT and hence the approximation orders 
needed for an accurate description. 

Higher-order variations are calculated by applying eq. (|37|) recursively or, alternatively, expand- 
ing the total rotation in ([3Tj) up to the required order. 

The second-order variation of the stress field a with the rotation X = p l ' X v can thus be expressed 
as the following symmetric 2-tensor: 

da\ 2 = [[a,X],X}=aXX-2XaX + XXa = p^-^^ ;i (38) 
which represents the second-order approximation space with coordinates p v ^ — p v p^ in the operator 

□2 

basis X vll = X y X^ = ggZgg^ , the action of which on a being self-apparent from this expression. 

The third- order variation of the stress field a with the rotation X = p l ' X v can be expressed as 
the following symmetric 2-tensor: 



S 3 a 

daL = \[[a,X},X},X} = aXXX - ZXaXX + SXXaX - XXXa = p v ^ - - - (39) 

This equation represents the third-order approximation space with coordinates p v ^ — p v p^p u in 

the operator basis X vfiu — X^X^X^ = ggugg^ggui ; the action of which on a is also straightforward. 

Higher order approximations are computed by following the same procedure. As a result one 

obtains a symmetric n-tensor of coordinates p u - K = p v ■ ... • p K in the basis X U ^ K = X v ■ . . . ■ X K = 

n ( 3 + n 1 \ 

g„ v d g„ K for each n-order. Also, each n-order variation encompasses a total of I I terms, 



n 

each of which can be factored into sets of n generators X v , where the resulting factors represent 



14 



the basis J„... K in the n-order approximation space. Therefore, each n-order is represented by the 
corresponding coordinates p v — K in the basis X V .„ K , and the coordinates represent a symmetric n- 

tensor with f ^ n ^ ~\ independent components. Note that the n-tensors of coordinates depend 

on the internal variables, p v = • • • , p v " K = p IA " K (£); tensors of coordinates can be properly 

taken as the position internal variables that describe anisotropy. 

Each variation has an equivalent interpretation, whatever its order, since transforming the axes 
of the reference system is equivalent to transforming the stresses involved — with a negative sign 
(sec Section 4) by virtue of the antisymmetric nature of the Lie bracket. 

The infinitesimal generators of conform to the following commutation relations: 

[X ai X b ]=s abc X c (40) 
where e a bc is the Levi-Civita symbol. These generators can be expressed as 2-tensors of the form 

\ 
0-1 
10/ 
1 \ 

(41) 
-10 0/ 
-1 \ 
10 
/ 

In summary, using the three generators of the group SOjl) in eq. (j4Tj) in conjunction with the 
Cauchy tensor a allows us to compute the variation of any order of the Cauchy tensor in order to 
obtain the different approximation to TSCT. 




8 n-Order variations of the stress tensor field 

This section specifically examines the variations of different order described in the previous section, 
which provide different approximations to TSCT. 

The first-order variation (|37|) of a with Xcan be expressed as 



[a,X]=pl 



P cr x 




2a xz 




®xx & 'zz 


-Vyz 





@xy 


: ~~ & 'zz 


®xy 


2&xz 



+ 



- P &yz + P (°~yy - °~xx 
p 1 2a yz - p 3 2cr xy 

P 1 (0- ZZ - O-yy) + P 2 a xy ~ p 3 O x 



) -P 1 0- X y+P 1 {0 xx 



a. 



P 1 {0-z 



-p 1 2a yz + p 2 2a x 



p\ 

p 2 o- xy - p 3 o 



(42) 
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which is a real symmetric traceless 2-tensor parameterized in the coordinates p v . The first-order 
variation for the particular case of the Cauchy diagonal tensor is 

/ p 3 (er 2 -cri) p 2 0i-cr 3 ) 

[cr } X]= p s {a 2 -ai) p 1 {a 3 - a 2 ) j (43) 

\p 2 (ai-a 3 ) p 1 ((T3-cr 2 ) 

and the second-order variation (|3"51) of a with X is 



[[ct, Jf],X] n = p 12 2a xy + p 13 2a xz - p 23 4a yz + p 33 2(a yy - a xx ) + p 22 2{a zz - <r xx ); 
[[a,X},X] 12 = ~{p n + p 22 + Ap 33 )a xy + p 23 3a xz + p 13 3a yz + p 12 (-2a zz + a yy + a xx ); 
[[a,X],X] 13 = p 23 3a xy - {p 11 + p 33 + Ap 22 )a xz + p 12 3a yz + p 13 {-2a yy + a xx + a zz ); 
[[a,X],X] 22 = 2p 12 a xy - p 13 Aa xz + p 23 2a yz - p 33 2(a yy - a xx ) + p lx 2{a zz - a yy ); [ > 

[[a,X],X] 23 = p 13 3a xy + p 12 3a xz + p 23 {-2a xx + a yy + a zz ) - (4p n + p 22 + p 33 )a yz ; 
[[a,X],X] 33 = -p 12 4a xy + p 13 2a xz + p 23 2a yz - p 22 2(a zz - a xx ) - p Y1 2(o zz - a yy ) 

which is a real symmetric traceless 2-tensor parameterized in the coordinates p u ^ — p v p^ . The 
specific expression for the Cauchy diagonal tensor is 



[[<r,X],X] 

p 33 2(a 2 -<7i) +p 22 2(a 3 -ay) p 12 {-2a 3 + a 2 + ai)) p 13 {~2a 2 + a x + a 3 ) 

p 12 (-2a 3 + a 2 +a 1 ) -p 33 2{a 2 -a 1 ) + p 11 2{a 3 ~a 2 ) p 23 (-2ai + a 2 + a 3 ) 

p 13 (-2a 2 + a 1 +a 3 ) p 23 {-2a x +a 2 + a 3 ) -p 22 2(er 3 - a{) - p n 2(er 3 - a 2 ) 



The third-order variation (|39f of a with X for the particular case of the Cauchy diagonal tensor 
is given by: 



a,X],X],X] u 


= p 123 6(a 2 - 


0-3); 








a,X],X],X]U 


= p 333 4(cri - 


cr 2 ) + 


p 223 (4 ( 7 1 - 3 ( T3- 


CT 2 ) +p 133 (-4cr 2 + 3cr 3 




a,X],X],X} 13 


= p 233 (-4a 1 


+ 3(T 2 


+ cr 3 ) +p 112 (4cr 3 


-3a 2 -a 1 )+p 222 A{a 3 


- 01); 


a,X],X],X} 22 


= P 123 6(a 3 - 










a,X],X],X} 23 


= p 133 (4a 2 - 


3(Ti — 


<7 3 ) +p m 4((T 2 - 


03) +p 122 (-4cr 3 + 3f7i 


+ ct 2 ); 


a,X],X],X} 33 


= p 123 6(<n - 


02) 









(46) 



which is a real symmetric traceless 2-tensor parameterized in the coordinates p v ^ = p v p^p u . 

The following series expansion encompasses those for all variation orders and defines the overall 
TSCT: 



da = da\ 1 + — da\ 2 + — da\ 3 + ... = [a,X] + — [[a,X],X] + — [[[a,X],X],X] + (47) 

Therefore, TSCT is represented by a symmetric traceless 2-tensor (j4"T|) that describes hydrostatic- 
pressure independent plasticity. If the solid concerned is completely isotropic, using the first term in 
(|47|) will suffice; if it is orthotropic, the first two will be enough. The plastic processes independent 
of the hydrostatic pressure are completely characterized by the stress dependences and tensors of 
coordinates p v , . . . 1 p v ' K — in their corresponding bases — in the terms for the different variation 
orders in fT?]) . 
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9 The norm space in plastic states 



The geometric study of plasticity conducted here involves examining hydrostatic-pressure indepen- 
dent plasticity in the tangent space to the Cauchy tensor (TSCT). As shown above, such a space 
can be represented by the scries expansion (|47[) , which has the typical form for a 2-tensor and gives 
both the spatial coordinates (degree of anisotropy) in a specific basis and the stress-dependence of 
hydrostatic-pressure independent plasticity. However, determining the plasticity potential (plasticity 
interaction) in this situation entails defining a specific metric for the 2-tensor (|47[) . To this extend, 
expression (l4"T|) — or the corresponding partial summation — is used to calculate the scalar norm g. 
The parameters introduced by the definition of the norm are restricted by the constraint that the 
plasticity function should be convex (third postulate) . 

Plasticity depends both on the externally applied stresses, a kl , and on the intrinsic properties of 
the material, p v , . . . , / o I/ " K , in a specific reference system or basis. Therefore, the scalar function g 
defining the norm will depend on both [i.e., g = g(a kl ; p v ' , p"'" K )]. The function g should provide 
the correlation between multiaxial stresses and hence indicate how the different components of the 
Cauchy tensor must be related. If (a kl )p is defined as the 2-tensor describing the plastic limit of 
the solid (stress internal variables), then such a limit will be reached when 

g(a kl ; (T, p"-») = g((a kl ) P ; pT, ....,>" ) (48) 

Usually, the tensor (o~ kl )p encompasses a single non-zero uniaxial component and thus coincides 
with op (the scalar uniaxial plastic limit that it is an internal variable). This allows the plasticity 
function to be rewritten as follows: 

f(a kl , (a kl ) P ;p\ ?"■«) = g(a kl ;p v , ..../-""-) - g((a kl ) P ; p" ..../-""-) (49) 

The scalar function g can be defined in terms of a matrix norm ||-|| such as L m (i.e. the usual 
m-norm for vectors); therefore, 

\\do-\r = (e?=i Ej 3 =i fan r) (m) 

g=\\da\r^f=\\da\r-\\(da) p \\ m 

where the constraints on the norm parameter (m S Sft, m > 1) ensure convexity in the plasticity 
yield criterion (third postulate) . The specific value of m depends on the properties of the particular 
material. Norm (|50p can be induced from the scalar product of second-order tensors, which is 
constructed from the following: 



< da, da' >= da kl da kl ' R T daR, R T da'R >= (Rlda kl Rf)(R k da' ' M R l 3 ) = da kl da kl ' (51) 

This definition, where rotations R\ and Rj are orthogonal (R T R — I), is invariant under rotations 
(R). The scalar product and its induced norm [eqs (|50|) and ([51]) ] afford measurement of angles and 
lengths in TSCT l[i7|). 

The definition must be completed with a scalar product of the coordinate tensors at the different 
independent orders (p v , . . . , p v " K ) of the form 

p»- K p x ....S = ( P n 2 ----(p K ) 2 5 v x ■■■■■■ (52) 
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The norm induced from this scalar product is L m for the corresponding n-coordinate tensor. 

In summary, plasticity yield criteria can be defined by the norms L m induced from the scalar 
products of stresses (I5T1) and coordinates (|52p in the stress and coordinate (position) spaces defined 
by the series expansion (|4"T)) . As a result of this construction, the space (|47p is normed. 



10 Analysis of plasticity yield criteria. Convexity 

The norms defined in Section 9 allow one to assess measurements of variations of any order with a 
view to establishing the plasticity yield criterion for each material in terms of its particular properties. 

Isotropic criteria can be analysed by using the norm L m for the first-order variation (|42|) . which 
gives 
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(53) 

With complete spatial isotropy (p 1 = p 2 = p 3 ) and invariance (rn = 2), the outcome is the Huber- 
Von Mises criterion (|13p . In the particular case of using independent plane stresses xy,yz,xz in 
each coordinate directions (0, 0, p 3 ), (p 1 , 0, 0), (0, p 2 , 0) respectively, taking (m = 2) and factoring in 
such a way as to make the criterion smooth (i.e., differentiable) leads to the Tresca criterion (|12p . 
With complete spatial isotropy (p 1 = p 2 — p 3 ), using the Cauchy tensor in principal stress (Cauchy 
diagonal tensor) gives the Hosford criterion (|17p . These are the most widely used basic plasticity 
yield criteria for materials with completely isotropic symmetry. 

Anisotropic criteria, which are applicable to orthotropic materials, can be analysed by using the 
norm L m for the second-order variation (|44p . which gives 
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(54) 



For the particular case of the Cauchy diagonal tensor, this leads to Hill's second criterion (flT 
Diagonalizing the coordinate tensor p v ^ and assuming (m = 2), we are led to Hill's first criterion 
(fl4|) . On the other hand, if both the coordinate tensor, p v>1 , and the Cauchy tensor are diagonalized, 
then one obtains the Logan-Hosford criterion l|18p . These criteria hold for anisotropic materials with 
orthotropic symmetry (i.e., orthotropic materials), a general description of which requires six spatial 
parameters (the six parameters of the symmetry tensor p 1 ^). 

The criteria of Barlat et al. [2] and Karafillis & Boyce [22] were originally established by 
transforming the stress tensor; this involved considering symmetry and anisotropy in the material. 
Within the framework of the proposed theory, the transformations required to describe these ef- 
fects in materials can be analysed via transformations of the coordinate tensors p v , . . . , p v - K from 
their starting bases (X v , . . . ,X„.... K ), to those needed for an accurate description of the material 
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concerned. Obviously, anisotropy makes solid materials basis-dependent (i.e., reference system- 
dependent); therefore, it grows with increasing order of the bases used. In completely isotropic 
materials, the plasticity yield criterion depends on the coordinates of a single scalar (p 1 — p 2 — p 3 ), 
i.e., it is absolutely independent of the particular coordinate system. 

In should be noted that Hill's second criterion, which provides an accurate description of the 
anomalous behaviour of some materials such as aluminium [18, 37], arises naturally within the frame- 
work of the proposed plasticity theory. Additionally, the proposed theory includes the anisotropy 
criterion of Hosford [25] that can explain the behaviour of fee (for m = 8) and bee (for m = 6) 
materials. 

With the norms defined above, equation (|47|) describes the most salient plasticity yield cri- 
teria reported so far as particular cases. Therefore, (I47p provides a self-contained description of 
hydrostatic-pressure independent plasticity (i.e., plasticity compliant with the proposed postulates 
and approximations). A higher-order analysis allows one to establish a generic criterion in terms of 
the following expanded series: 
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(55) 

where i ^ j ^ k ^ i and the indices can range over x, y, z; f is the permutation index. Each term in 
(|55[) has an associated coefficient u n i f,v,w,k that can be explicitly computed within the framework 
of the proposed unified theory, but has been excluded in the expanded from (|55[) for simplicity. Note 
that positive integer parameter n = n' + 1 in eq. f|55j) represents the n-order approximation. 

Stress convexity in (|55|) can be assessed by assuming that a linear combination of convex functions 
with positive coefficients is a convex function itself [36]. The coefficients u n 'f,v,w,k are in fact all 
positive. Therefore, it will suffice to show that the stress-dependence of each individual coefficient 



in (|55|) is convex (i.e., that all stresses in the terms 



IK 



CTfcfc 



are convex) . Convexity in a function T of three variables can be expressed 



T(cia u + c 2 o-jj + c 3 o k k) < ciT(o- ti ) + C2T(ajj) + c 3 T(a k k) 



(56) 



where ci, 02,03 € 5? are positive or zero and the following condition holds: c± + C2 + C3 — 1. The 



convexity of \a xy 



\o~ xz \ m can be immediately confirmed by applying (|56p to the particular 



case of a single variable. That of 



2" cr, 



can also be easily inferred from the 



(2» - l)a n - a kk 

triangular inequality (see Appendix) fulfilling this L m norm as defined in the stress vector space 
(o~ xx , Cyy, a zz ). Note that, for L m to actually be the norm, m £ 5ft, m > 1, which is the constraint to 
be imposed on (|55p for stresses to be convex. Under these conditions, the plasticity yield criterion 
(|55p will fulfil the third postulate. Within the framework of the proposed theory, the general criterion 
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(|55p can predict plastic behaviour in solids. Therefore, it can be of use in future experimental and 
theory studies on this topic. 

11 Hydrostatic pressure-dependent plasticity yield criteria 

The hydrostatic pressure-independent definitions obtained in the previous sections invariably include 
traceless 2-tensors. In fact, this mathematical constraint ensures that the ensuing criteria will be 
hydrostatic-pressure independent. As noted earlier, the normal space to the Cauchy tensor describes 
the hydrostatic pressure-dependence of plasticity. As a result, considering the effects of hydrostatic 
pressure in plasticity yield criteria entails including a c(tTA/) fc = c (cr xx + a yy + cr zz ) b dependent 
function to describe the hydrostatic behaviour of the material with specific coefficients b, c [10, 11, 
7, 9]. A more detailed study of the hydrostatic pressure-dependence in the framework of unified 
theory will be the subject of future work. 

12 Conclusions 

In this paper we propose a new approach to the theory of macroscopic plasticity, based on a geometri- 
cal ansatz, that unifies the different plasticity yield criteria. The method is based on decomposing the 
Cauchy tensor into its component stress spaces: tangent and normal spaces. Hydrostatic-pressure 
independent yield criteria are located on the tangent space to the Cauchy tensor (TSCT), whereas 
hydrostatic-pressure dependent yield criteria are located on the tangent and normal spaces to the 
Cauchy tensor (TSCT and NSCT). The analysis of the corresponding spaces requires the use of Lie 
groups of transformations. 

This work focuses on the tangent space to Cauchy tensor (TSCT) and analyzes the hydrostatic- 
pressure independent yield criteria as orthogonal transformations. The measure of this tangent 
space is taken as a mathematical norm that guarantees the convexity of the yield criteria. From 
this measure we obtained the convex-stress dependence for both isotropic yield criteria (Tresca 
1864 and Von Mises 1913) and anisotropic yield criteria (Hill 1948 and Hill 1979). In addition, 
a stress series expansion is obtained that allows a detailed study of the mechanical properties in 
materials: anisotropy, hardcning-softcning, etc. The tensor coefficients resulting from orthogonal 
transformations are dependent on internal variables. These tensor coefficients allow a deep analysis 
of mechanical properties in materials. This unified theory simplifies the great existing amount of 
yield criteria and gives a physical meaning to the anisotropy coefficients. 

This unified theory also allows an analytical study of anisotropy (dependency of /on £) and 
hardening-softening (dependency of / on £) by means of a set of tensor coefficients. These coefficients 
are the directions of macroscopic plasticity (by points) in materials. A connection is suggested 
between macroscopic transformations (unified theory) and the microscopic processes (dislocation 
theory). To this extent, the concept of macroscopic point, that constitutes a statistical measure of 
microscopic processes, is introduced. The concept of the point connects two scales with different 
physical laws. 

This unified approach has a series of advantages. At first, it proposes the stress dependencies 
of yield criteria, starting from very simple postulates and with a unified vision of the underlying 
phenomena. Further, it allows to know the macroscopic plasticity directions in the solid points, 
by tensor coefficients of orthogonal transformations. It also provides a manner to analyze the 
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anisotropy and hardcning-softening in materials. Finally, it enables us to establish a connection 
between microscopic plasticity and macroscopic plasticity. 

In future work it is important to extend the unified theory to treat the very important question 
about the hydrostatic pressure-dependent plasticity (important in composites and granular materi- 
als); for this it is necessary to analyse dilatation transformations. In this context, it is of fundamental 
importance to make a deep study of hardening-softening and the time-evolution of plasticity in ma- 
terials; for the latter aspect it is necessary to determine an adequate time-transformation that gives 
a correct parameterization of the solid time-evolution. 
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Notation 

(a) Indices 

1. Solid point position: A ranging over x, y, z 

2. Cauchy tensor: i,j, k, l,p, q ranging over x, y, z 

3. Rotation group: v, n, w, k, x, £ ranging over 1, 2, 3 

4. Levi-Civita tensor: a, 6, c ranging over 1, 2, 3 

5. Elastic-plastic limit: P 

6. Angle equivalence: (fl 1 ,*? 2 ,*? 3 ) = {a, (3, 7) 

7. Cauchy tensor equivalence: a = a kl (contravariant); a = C/t; (covariant) 

8. Principal stresses of the Cauchy tensor: a\, <J2, &3 

9. Hydrostatic pressure: gm 
10. Internal variables: £ = £ Q 

(b) Tensors 

The Einstein summation convention is used throughout. Covariant tensor notation is systematically 
employed from Section 4 onwards. 
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Appendix: Convexity of general plasticity function 



In order to prove the convexity of the plasticity function (section 10), it is necessary to show the 



convexity of the terms 



2" a, 



(2" 



IK 



0~kk 



with m > 1, m G 5ft and n' a positive integer 



or zero. If the triangular inequality is taken we obtain 



2" <7,: 



(2"' - l) <r S 



0~kk 



< 2 



Dividing (A) by 2 m ("' +1 ) results in 



|o-«r + (2"' - l) \a J3 \ m + \a kk \ m (A) 



2" + 1 
Wkk 



'J J 



2 „' + i °"fcfc 



< 



1 

2"' + ! 



'Ml 



+ 



< 2 



(B) 



2 m(n'+l) 

Thus taking T(x) = \x\ m , the condition of the convexity is fulfilled 
T{ciau + c 2 (Tjj + c 3 a k k) < CiT(a u ) + c 2 T{a :jj ) + c 3 T(a kk ) (C), 



where c\ 



and Ci + c 2 + C3 = 1. So the convexity of the plasticity function has been proved. 



} +1 and n' is a positive integer or zero; with c\, c 2 , C3 > 
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Figure captions 



Figure 1. Reference systems for the overall solid {0,x x — (x, y, z)} and for each point p in it 
{O p ,x x — (x, y, z)}. Each point has an associated Cauchy tensor acting on it (see Figure 2). 

Figure 2. Representation of the Cauchy tensor a kl and its rotational possibilities about the Eulcr 
angles Qe = &e{oi,(3, 7) in the reference system {O p ,x A }. The sphere surface in the figure is a 
particular making of the S'Ofjnj) rotation group. 

Figure 3. Action of the .SO© rotation group on the Cauchy tensor a kl for a partial rotation f3 
about the y-axis, represented by the reference system transformation {O p ,x A } — > {O p ,x A }. 

Figure 4. Arbitrary representation of TSCT as an ellipsoid. The tangent space TSO^ pl _, that 
provides a first-order approximation of TSCT to a reduced neighbourhood of is shown. 

Figure 5. Arbitrary representation of TSCT and its approximations of different order to an in- 
finitesimal neighbourhood. The higher orders lead to local infinitesimal neighbourhood that provide 
a more accurate description of the overall structure of TSCT. 
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Figure 1: Reference systems for the overall solid {0,x A } = ({x, y, z})} and for each point p in it 
{O p ,x x } — ({x, y, z})}. Each point has an associated Cauchy tensor acting on it (see Figure 2). 
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Figure 2: Representation of the Cauchy tensor a and its rotational possibilities about the Euler 
angles Qe = 0e(o!,/3, 7) in the reference system {O p ,x A }. The sphere surface in the figure is a 
particular making of the SO(1?) rotation group. 
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c/dy 




Figure 4: Arbitrary representation of TSCT as an ellipsoid. The tangent space TSO^ pl + that 
provides a first-order approximation of TSCT to a reduced neighbourhood of is shown. 
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Figure 5: Arbitrary representation of TSCT and its approximations of different order to an in- 
finitesimal neighbourhood. The higher orders lead to local infinitesimal neighbourhood that provide 
a more accurate description of the overall structure of TSCT. 
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